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ABSTRACT 

Misaligned accretion disks surrounding rotating compact objects experience a torque due to the Lense- 
Thirring effect, which leads to precession of the inner disk. It has been suggested that this effect could 
be responsible for some low frequency Quasi-Periodic Oscillations observed in the X-ray lightcurves of 
neutron star and galactic black hole systems. We investigate this possibility via time-dependent calcu- 
lations of the response of the inner disk to impulsive perturbations for both Newtonian point mass and 
Paczynski-Wiita potentials, and compare the results to the predictions of the linearized twisted accretion 
disk equations. For most of a wide range of disk models that we have considered, the combination of 
differential precession and viscosity causes the warps to decay extremely rapidly. Moreover, at least for 
relatively slowly rotating objects, linear calculations in a Newtonian point mass potential provide a good 
measure of the damping rate, provided only that the timescale for precession is much shorter than the 
viscous time in the inner disk. The typically rapid decay rates suggest that coherent precession of a 
fluid disk would not be observable, though it remains possible that the damping rate of warp in the disk 
could be low enough to permit weakly coherent signals from Lense-Thirring precession. 

Subject headings: accretion, accretion disks — hydrodynamics — black hole physics - relativity ~ stars: 
neutron - X-rays: general 

Animation of warp evolution is available at: http://www.cita.utoronto.ca/~arniitage/lense_thirriiig.litml 



1. INTRODUCTION 

Tilted orbits about a rotating object experience a torque 
due to the general relativistic Lense-Thirring effect (Lense 
& Thirring 1918), which causes the plane of the orbit to 
precess. For a fluid accretion disk, the differential pre- 
cession with radius causes stresses and dissipation. If the 
torque is strong enough compared to the internal viscous 
forces, the result is that the inner regions of the disk are 
forced to align with the spin of the central neutron star or 
black hole (Bardeen & Petterson 1975). 

This process has several important effects. If the an- 
gular momentum of disk material at large radius is mis- 
aligned with respect to the rotation axis of the central ob- 
ject, then alignment of the inner regions by the Bardeen- 
Petterson effect implies that the disk must possess a large- 
scale warped or twisted shape. It is then obvious that this 
will both modify the emergent spectrum, and determine 
the direction of jets accelerated from the inner disk re- 
gion. The consequences of this for disks and jets in Ac- 
tive Galactic Nuclei are well known, and have been dis- 
cussed for some time (Bardeen & Petterson 1975; Rees 
1978; Scheuer & Feiler 1996; Natarajan & Pringle 1998; 
Natarajan & Armitage 1999). 

More recently, it has been suggested that Lense-Thirring 
precession may be directly observable in the form of Quasi- 
Periodic Oscillations (QPOs) in the X-ray lightcurves of 
low- mass X-ray binaries (Stella & Vietri 1998), including 
galactic black hole candidates (Cui, Zhang & Chen 1998). 
If confirmed, this identification promises to provide con- 
straints on the equation of state of material at nuclear 
densities, and on the spin parameter of stellar mass black 
holes. It would also constitute one of the few arenas where 



the Lense-Thirring effect might be detectable - the mea- 
surement of this precession for orbits around the Earth 
(Ciufolini et al. 1998) is, itself, controversial. However, 
other explanations for the observed QPOs remain viable. 
For neutron stars especially, the possibilities for interac- 
tion between the disk and a magnetosphere are legion, 
and in general remain poorly understood (see e.g. Miller 
& Stone 1997; Spruit & Taam 1993; Ghosh & Lamb 1979). 
There are also other relativistic effects possible for purely 
flat, fluid disks around black holes, which can produce os- 
cillations with frequencies in the right range (e.g. Nowak 
et al. 1997). 

Estimates of the frequency of Lense-Thirring precession 
for parameters appropriate to neutron stars are, at least to 
an order of magnitude, in agreement with the observed fre- 
quencies of some QPOs (Stella & Vietri 1998). The main 
theoretical uncertainty in applying the model is then the 
damping rate of perturbations excited in the inner disk. 
This was computed by Markovic & Lamb (1998), who 
considered the global modes of the disk described by the 
linearized version of the twisted disk evolution equations. 
They found that the combination of differential precession 
and disk viscosity was generally highly deleterious to the 
survival of warps in the disk. The damping rates they 
obtained were sufficiently rapid as to cast serious doubt 
on the Lense-Thirring interpretation of QPOs, although 
some relatively weakly damped modes were found when 
the torques from radiation were also included. In response, 
Vietri & Stella (1998) suggested that damping of verti- 
cal motions might be much weaker than that assumed by 
Markovic & Lamb (1998), due either to a radically modi- 
fied disk structure, or a supposed more general suppression 
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of damping in the inner, aligned portions of the disk. 

In this paper, we compute numerically the response of 
a viscous disk to finite amplitude perturbations in the in- 
nermost regions where the Lense-Thirring torque is strong. 
Following Markovic & Lamb (1998), we employ the equa- 
tions derived by Papaloizou & Pringle (1983) that describe 
the hydrodynamic response of a thin, viscous disk that 
evolves essentially diffusively. Our approach is to com- 
mence with a steady, flat disk, perturb the inner disk 
impulsively to generate a warp, and then evolve the full 
twisted disk equations to follow the warp as it precesses 
and decays. We compare our results both with the linear 
evolution investigated previously, and with the evolution 
in a Paczynski-Wiita (1980) potential that models some of 
the non-Newtonian effects expected close to neutron stars 
and black holes. 

The outline of the paper is as follows. We set out the 
equations solved in §2, and estimate the parameters ap- 
propriate for accretion disks in X-ray binaries. In §3 we 
present results from a grid of models with varying viscosity 
laws, potentials, and approximations to the twist equation. 
§4 summarizes our conclusions. 

2. EVOLUTION EQUATIONS FOR A TWISTED DISK 

The evolution of a geometrically thin, warped accre- 
tion disk in the diffusive regime can be described using 
the equations developed by Papaloizou & Pringle (1983). 
In terms of the angular momentum density of the disk 
L(i?, i), the governing equation is (Pringle 1992), 
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Here L = [GMf/'^Y.R^/'^i, S is the surface density, and I 
is a unit vector normal to the disk surface. The angular 
velocity in the disk at radius R is ^{R), and fi = dO/di?. 
We consider Keplerian orbits in two possible forms for the 
potential, a Newtonian point mass potential. 
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and a pseudo-Newtonian potential (Paczynski & Wiita 
1980), 

^--^ (3) 

R — Rg 

where Rg = 2GM/c^. The latter potential reproduces cor- 
rectly some features of orbits around compact objects - for 
example the existence and location of an innermost stable 
orbit - though we include it here mainly to gauge the sen- 
sitivity of equation (1) to changes in the angular velocity 
profile, vi and 1^2 are "viscosities" corresponding to the az- 
imuthal and vertical shear respectively, and the last term 
on the right hand side of the equation is the Lense-Thirring 
torque. For a black hole with angular momentum J, for 



example, flp = uJp/R^ = 2JG/c^i?^. We note that this 
fundamentally Newtonian approach is appropriate only for 
black holes with modest values of the spin parameter. In 
particular, the behavior of disks around rapidly rotating, 
Kerr black holes, is formally not addressed by these cal- 
culations. Although there are no direct measurements of 
black hole spins, observations are suggestive of a broad 
range of spin parameters for black holes in X-ray binaries 
(Zhang, Cui & Chen 1997). 

Equation (1) represents a vast simplification to the 
general equations for the time dependent evolution of a 
warped accretion disk, which in general require a full hy- 
drodynamic treatment (see e.g. the simulations of Lar- 
wood et al. 1996). The conditions under which equation 
(1) is valid have been investigated in detail (Papaloizou & 
Pringle 1983; Papaloizou & Lin 1994; Ogilvie 1999), with 
the result that diffusive evolution is likely to apply to thin 
disks in Active Galactic Nuclei and, probably. X-ray bi- 
naries (Pringle 1999), which are the systems of interest 
here. Conversely wave-like evolution is likely to dominate 
for disks around young stars. 

Equation (1) can be derived either from consideration of 
mass and angular momentum conservation of two neigh- 
boring annuli in the disk sharing their angular momentum 
via viscosity (Papaloizou & Pringle 1983; Pringle 1992), 
or via a hydrodynamic treatment of the warped disk valid 
in the linear regime (Papaloizou & Pringle 1983). More 
recent analysis (Ogilvie 1999) has shown that an equa- 
tion of this type is valid even for strongly warped disks, 
though, in general, neither vi nor 1^2 so defined are equal 
to the usual kinematic viscosity v of the Navier-Stokes 
equation and planar disk theory. In fact, the relation be- 
tween these quantities for an arbitrary warp is extremely 
complex (Ogilvie 1999). For this paper, we will take the 
ratio of vi to V2 to be a free parameter, and examine the 
influence of this choice on the decay of warps in the inner 
disk. Clearly, models with low 1^2 viscosity are the most 
favorable for sustaining observable precession of the inner 
disk regions. 

There are two relevant timcscales in this problem, the 
precession time and the viscous timescale corresponding 
to the vertical shear, which we define as, 
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The balance between these timescales determines whether 
the Bardeen-Petterson effect is able to align the inner disk 
with the spin of the central object. Roughly speaking, this 
will occur if tprec <, tu2 at R = Rin (Kumar & Pringle 1985; 
Scheuer & Feiler 1996; Natarajan & Armitage 1999). 

This ratio of timescales, t-pT:cc/tv2 1 can be estimated as 
follows. For definiteness, we assume that the disk sur- 
rounds a black hole with spin parameter a, for which 
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and the corresponding precession timescale is. 
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If the J^i disk viscosity is parameterized via the Shakura- 
Sunyaev (1973) prescription, i^i = ac^/fi, where a is a 
dimensionless parameter and Cs is the local sound speed, 
then ti,^ ~ B? /vi is given by. 
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where H is the disk scale height. We then have, 
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where we have adopted rough estimates for a and H/R. 
Obviously, there are large variations possible in all of these 
parameters, for example reducing H/R to 10^^ would de- 
crease tp-cec/tvi by two orders of magnitude. However, if 
t^2 ~ J^i, then this simple estimate suggests that the pro- 
cessional timescale in the inner disk is likely to be perhaps 
two or three orders of magnitude shorter than the local 
viscous timescale. Accordingly, we choose parameters for 
our calculations that match this regime, and the a priori 
more favorable one (for observing disk precession) where 
the viscosity, and hence the damping of warps, is weaker. 

We solve equation (1) using the explicit finite differ- 
ence method described by Pringle (1992), with zero torque 
boundary conditions at i?in and i?out- The outer bound- 
ary condition is applied at a radius large enough as to 
not affect the inner disk evolution over the time period of 
interest. Typically, i?out/-Rin = 16 suffices. Moderately 
high numerical resolution is needed to resolve the strongly 
twisted disk configurations that develop, we use 400 to 800 
grid points logarithmically spaced between _/?!„ and i?out- 
Tests show that this resolution is more than adequate. 

For our calculations we adopt units in which i?i„ = 1, 
and |fij,(i?in)| = 8. We take the viscosities to be power 
laws in radius, vi — i^ioR , V2 ~ 1^20 R, and consider 
cases where 1^20/^10 ~ 5,1,0.2. We consider an initially 
aligned, planar disk, apply a simple warp perturbation to 
the inclination. 
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and follow the perturbed disk as the warp decays and pro- 
cesses. We note that this form of perturbation, in which 
there is initially no twist in the disk warp, generally decays 
slower than a twisted configuration for which the radial 
scale on which components of I change is smaller. 

These initial conditions are appropriate if the outer ac- 
cretion disk is aligned with the equatorial plane of the 
spinning black hole. Alternatively, the outer disk could 
be misaligned, either as a consequence of warping insta- 
bilities (Pringle 1996; Maloncy & Begelman 1997; Schandl 
& Meyer 1994), or because the black hole retains an ini- 
tially misaligned angular momentum vector (King & Kolb 
1999). However, even in these cases, we expect that the 
inner disk will be aligned with the spin of the hole by the 
Bardeen-Petterson (1975) effect. For standard models of 
the disk viscosity, this alignment radius is large compared 
to the region of the disk from which QPOs are expected 
to originate. This conclusion is reasonably robust, at least 
for steady disks, since it follows simply from the relative 



timescale for precession versus viscous evolution. Further- 
more, we show later that models in which alignment would 
not occur (those with relatively large values of the disk 
viscosity), are those in which perturbations damp most 
rapidly. 

3. RESULTS 
3.1. Newtonian point mass potential 

Fig. 1 shows the decay of a warp in the inner regions of 
the disk. The parameters of the initial perturbation were: 
Rp = 2, ARp = 0.2, and Azq = 0.1. The azimuthal an- 
gle of the ascending node, 7, was initially constant with 
radius, i.e. the initial state had no twist. The time slices 
are plotted at intervals of At = 5, which is a little less 
than the precession time, tpicc = 27r, at Rp. The disk vis- 
cosity was taken to be i^i = 1^2 = lO^^i?'^/^, which gives 
tp,J/t^^ = 4 X 10^^ 

From the Figure, it can be seen that the warp decays 
rapidly, with the peak inclination decaying by roughly an 
order of magnitude per precession time. As it does so, the 
warp diffuses radially, and the differential precession with 
radius induces a strong twist in the disk shape. This is 
shown in Fig. 2, which plots the shape of the disk at the 
instant when the perturbation is applied and after a time 
Ai = 15, around two precession timescales later. The warp 
becomes increasingly twisted with time, which reduces the 
characteristic radial lengthscales for variations in Z, and 
thereby shortens the timescale (~ AR^ /V2) required for 
diffusion to flatten the disk. 

Fig. 3 shows the temporal behavior of the warp for a 
grid of disk models with varying assumed forms for the 
viscosity. We consider viscosity laws with 5 = 3/2 and 
(5 = 0, for a range of 1^10 and 1^20- 

For the cases where the viscosity is weak (i.e. low 1^20), 
the decay of the warp depends almost exclusively on the 
value of 1^2 for the range of i^i considered here. The ini- 
tial rate of decay is found to be proportional to 1^2^, but 
this phase lasts for at most a few precession times. Sub- 
sequently, the disk inclination is found to decline more 
rapidly with time, so that at the end of the runs the rate 
of decay d\ogi/dt is similar for the various choices of 1^2- 
The runs with differing S are very similar, the main distinc- 
tion being the trivial one that the i5 = calculations have 
lower viscosity at i? = _Rp due to our choice of normaliza- 
tion. Thus, the decay rates are correspondingly lower. 

Our strongest viscosity runs show some dependence of 
the damping rate on i^i as well as 1^2- The model with 
large v\/v2 exhibits substantially faster decay. This is due 
to the larger radial velocity advecting the perturbation in- 
wards, which causes a more rapid decline in the inclination 
at the fixed radius Rp. 

3.2. Precession rate 

Fig. 4 shows the precession rate, d7/dt, for the annu- 
lus at i? = 2 in the S = 3/2 models. The Lense-Thirring 
precession rate at this radius for our choice of units corre- 
sponds to (dj/dt) = 1, and this is the initial rate of pre- 
cession for the warped disk in all the models. Thereafter, 
the precession rate declines as the increasingly twisted disk 
generates internal viscous torques that oppose further dis- 
tortion of the disk shape, with the decay rate of the pre- 
cession varying with 1^2 in the same manner as the disk 
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inclination. Comparing Figs. 3 and 4, we find that the 
precession of the disk annulus remains within 10% of its 
initial value up to the epoch when the disk inclination has 
decayed by roughly an order of magnitude. For all the 
models considered here this is a short timespan, of the 
order of a few precession periods. 

Converting these results to a prediction of the coherence 
(or lack thereof) of a Lense-Thirring feature in the power 
spectrum of X-ray binaries would require detailed knowl- 
edge of the emission mechanisms producing the X-ray flux, 
which is generally lacking. However, it is clear that for the 
models considered here, the perturbed disk is able to main- 
tain a coherent precessional motion, at a fixed radius, for 
at most a few precession periods. This might be sufficient 
to generate weakly coherent oscillations in the disk emis- 
sivity (and some of the observed QPOs are indeed of this 
type), but would be unable to lead to oscillations with a 
high quality factor. Moreover, this estimate most proba- 
bly overstates the coherence of the X-ray signal, since any 
plausible emission mechanism will sample a range of radii 
which will have different precession frequencies. 

3.3. Comparison with linearized equations 

It is of interest to compare our results with the compre- 
hensive analysis of the linearized twist equations, includ- 
ing Lense-Thirring precession, presented by Markovic & 
Lamb (1998). For small warp amplitudes, where /3 ^ 1 
and Rf3 <C 1 (here (3 is the local angle of tilt), the surface 
density of the disk remains a fixed function of radius, and 
the warp can be completely described by two components 
of the local tilt vector, l^ and ly. This simplifies equa- 
tion (1) significantly, and the resulting equation can be 
expressed in terms of a single complex quantity W = j3e^^ . 
For V2 oc B?/^ the evolution of W is described by. 
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Markovic & Lamb (1998) proceeded to solve for the modes 
of this equation, and the corresponding equation including 
terms representing radiative torques. For our purposes of 
comparing results with the previous calculations, we in- 
stead solve equation (10) as an initial value problem with 
the initial conditions defined by equation (9). This can 
straightforwardly be done using identical methods to those 
employed for the full evolution equation. 

Figure 5 shows the comparison between the disk evolu- 
tion computed with equations (1) and (10). Two choices 
of viscosity are used, the 5 = 3/2, j^io = V20 = 10~^ model 
discussed already, and a model where the viscous time is an 
order of magnitude longer. Considering this second case 
first, for these parameters, where iproc ^ tv, the decay 
of the warp is essentially identical between the two cal- 
culations, both in the rate and the detailed warp shape. 
This regime is the only one where there is some possibil- 
ity of precession surviving for long enough to be poten- 
tially observable, and for warps of this (relatively modest) 
amplitude we therefore conclude that the linear approach 
provides an excellent approximation to the evolution de- 
scribed by equation (1). For the model where tprec and ti, 
are more comparable (the upper panel in Fig. 5), there are 
some differences between the detailed shape of the warp 
in the two calculations. These are due to changes in the 



surface density of the inner disk in response to the warp, 
which are ignored in the linear treatment. Even in this 
case though the decay rates are found to be very similar, 
and there is certainly no evidence that including the extra 
terms in equation (1) can cause significantly slower decay 
than is obtained in the linear regime. 

3.4. Paczynski-Wiita potential 

A rather analogous situation to the preceeding Section 
is found for the comparison between the Newtonian point 
mass and Paczynski-Wiita potentials. Figure 6 shows the 
evolution of disk inclination as a function of radius for the 
two potentials, using the same parameters as for Fig. 1. 
There are significant differences between the shape of the 
warp in the two calculations. At late times the warp in 
a Paczynski-Wiita potential sustains a higher amplitude 
at the inner disk edge than in the Newtonian case, and in 
principle this would be more favorable for producing an 
observable signal of disk precession. This is attributable 
to the lower surface density (and correspondingly more 
rapid radial inflow for a given disk viscosity) at small radii 
in the non-Newtonian potential. However, the damping 
rates for the two calculations differ only to a qualitatively 
insignificant degree - for this choice of parameters both 
warps decay extremely rapidly. Moreover, for less viscous 
disks, where the warps decay more slowly, the differences 
between the evolution in the two potentials quickly become 
much smaller. 

3.5. Radiation induced warping 

In the preceeding discussion, we have investigated the 
damping rates of warps in disks whose evolution is driven 
solely by induced precession and internal disk viscosity. 
The influence of radiation forces caused by re-emission of 
radiation from a central source of luminosity is known to 
lead to growing warping modes of the disk (Pringle 1996) 
in some circumstances, so it is of interest to consider how 
they may alter the evolution of the inner disk. The flrst 
point to note is that the large scale warping modes dis- 
cussed by Maloney, Begelman & Pringle (1996; see also 
Maloney & Begelman 1997; Maloney, Begelman & Nowak 
1998; Wijers & Pringle 1999) are probably irrelevant to the 
discussion. These modes grow and precess on timescales 
that are of the order of the viscous timescale of the outer 
edge of the disk, typically tens to hundreds of days for 
X-ray binaries. The most promising circumstances for ob- 
serving Lense-Thirring precession require tproc ■^ ^1^2 in 
the inner disk, in which case any outer warp will be flat- 
tened into the equatorial plane well outside the innermost 
regions. However, the low surface density of the inner 
regions, if they remain optically thick, allow for the pos- 
sibility that radiation reaction forces could play an addi- 
tional signiflcant role there. Indeed, consideration of the 
linearized twist equations, including the radiative torque 
(Markovic & Lamb 1998), suggests that the decay rate of 
some twisted disk modes could be substantially reduced 
over the case where radiation is neglected. 

Although highly suggestive, we feel that the influence 
of radiation reaction forces on the survival of inner disk 
precession is likely to be considerably more subtle, ow- 
ing primarily to the importance of shadowing. Shadow- 
ing of parts of the outer disk by twists at smaller radii 



Armitage & Natarajan 



is an inherently non-linear effect that provides additional 
coupling between different radii in the disk. For a disk 
illuminated by a point source at the origin, it modifies the 
growth rate even for negligible disk inclinations, where the 
linear approximation would otherwise be extremely accu- 
rate. We have found that the sense of the modification 
is typically to reduce the growth rate for disks that are 
subject to radiation warping at large radius, and we an- 
ticipate that there would be equally significant effects for 
the much more tightly wound modes present at small radii 
when the effects of Lense-Thirring precession are included. 
Additionally, at radii of just a few Rg, the assumption of a 
central illuminating point source (which is an excellent ap- 
proximation for the usual radiation driven modes at very 
large disk radii) breaks down, and the real radiation field 
is likely to be much more complex. It is unclear to what 
degree these effects alter the damping rates computed by 
Markovic & Lamb (1998), but clearly further investigation 
in this area is merited. 

4. DISCUSSION 

In this paper, we have presented results from a time- 
dependent treatment of the evolution of perturbations in 
the inner regions of accretion disks which are aligned with 
the spin axis of the accretor by the Bardeen-Petterson ef- 
fect. We find that such perturbations decay rapidly. For 
a weak viscosity (one for which the viscous timescale cor- 
responding to i^2 is much greater than the local precession 
timescale), the inner parts of the disk would be expected 
to be accurately aligned with the rotation axis. An impul- 
sively generated warp then decays at a rate that depends 
solely on the viscosity acting on the {R, z) stress, conven- 
tionally parameterized as 1^2 (Papaloizou & Pringle 1983). 
The initial damping rate is proportional to 1^2, but after 
a few precession periods the damping is greatly enhanced 
as a consequence of the strong differential precession and 
twisting of the disk shape. This twisting rapidly reduces 
the radial lengthscale over which disk viscosity must act 
in order to flatten the disk into the equatorial plane, and 
leads to rapid damping. 

Our calculations are based on a numerical treatment 
of the full twisted disk equations derived by Papaloizou 
& Pringle (1983), which can readily be extended to arbi- 
trary rotation laws (Pringle 1992). We have compared the 
evolution of warps in potentials arising from a Newtonian 
point mass, and a Paczynski-Wiita potential that mimics 
some features of hydrodynamics in a Schwarzschild met- 
ric. For disk parameters chosen such that warps are at 
least relatively long lived, we find that the decay rate of 
warps in these two potentials are very similar, and accu- 
rately described by a linearized treatment of the twisted 
disk equations. Disks where the inflow velocity is larger 
exhibit some qualitative differences, both between the two 
potentials and when compared to the linear results, but 
these disks show even faster damping than expected on 
the basis of the value of V2 alone. They are not favorable 



cases for producing long lived, observable, precession. Ad- 
vection dominated flows (e.g. Narayan & Yi 1994), where 
the radial inflow velocity is typically a significant fraction 
of the orbital velocity, would be extreme examples of this 
effect. Such disks seem even less likely to be able to sup- 
port coherent Lense-Thirring precession. 

These calculations, and those of Markovic & Lamb 
(1998), assume that the accreting material in the vicin- 
ity of the black hole can be treated as a fluid accretion 
disk. This would be wrong if the inner disk was prone 
to breaking into blobs, for example via the development 
of the Lightman-Eardley (1974) or other instabilities (e.g. 
Ghosh 1998). Several authors have suggested that the in- 
ner disk in X-ray binaries is, indeed, susceptible to such 
instabilities (e.g. Taam & Lin 1984; Lasota & Pelat 1991; 
Cannizzo 1996). The behavior of accreting gas in such a 
model could be very different from that considered here, 
and is beyond the scope of this paper. Further study of 
such instabilities, in the context of currently favored mag- 
netohydrodynamic models for the origin of the disk vis- 
cosity, would be valuable. 

The strength of the viscosity in planar accretion disks 
is constrained, albeit poorly, by a variety of observational 
methods (see e.g. Cannizzo, Chen & Livio 1995). By con- 
trast, there are essentially no constraints whatsoever on 
how efficiently a disk can damp warp, especially at the 
small radii that are of relevance to Lense-Thirring pre- 
cession (at large radii weak constraints are possible from 
observing which systems may be prone to radiation in- 
duced warping). This introduces a major uncertainty in 
any attempt to apply the results to interpreting the possi- 
ble origins of QPOs in X-ray binaries. However, if vi ~ 1^2, 
it seems clear that warps in the inner disk are damped on a 
timescale of the order of a single precession period, or less. 
This is the same result as was found by Markovic & Lamb 
(1998) from a linear analysis, and it suggests that preces- 
sion is unlikely to be the origin of low-frequency QPOs in 
disks around neutron stars and galactic black hole candi- 
date sources. The important caveat is that we know of 
no definite reason, cither from observational constraints 
or from theoretical arguments, why 1^2 might not be very 
much smaller that vi in the inner disk. In this regime, dif- 
ferential precession still leads to strong twisting and even- 
tual rapid decay of the warp. However, our lowest 1^2 runs 
maintain a reasonable warp amplitude for ^ 10 precession 
periods, along with a roughly constant precession rate. 
This could lead to at least weakly coherent signals that 
have as their origin Lense-Thirring precession. Better the- 
oretical understanding of the hydrodynamics of disk warps 
will be needed to investigate this possibility. Recent work 
(Ogilvie 1999) has made progress towards this goal. 

We thank Jim Pringle for making available his numeri- 
cal disk evolution code, Gordon Ogilvie for sharing results 
in advance of publication, and the referee Wei Cui for a 
prompt and very helpful report. 
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Fig. 1. — Decay of the warp for the run with ui = U2 = 10 ^B?''^ . From top downwards, time slices plot the inclination angle in radians 
at intervals of At = 5. The time units are such that |f2p|^^ = 1 at i? = 2. 
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Fig. 2. — Disk shape for the run with ui = U2 = W^^B?'^ . The left image shows the initial warp perturbation at t = 0, the right image 
the warp at t = 15. Note that the amplitude of the warp has been vastly exaggerated in the second image to show the shape of the rapidly 
decaying perturbation. 
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Fig. 3. — Damping of the disk inclination i at R = 2Hin for varying disk viscosity. Upper panel is for 5 = 3/2, lower panel is for <5 = 0. Solid 
curves are for vxjv^ = 1, dashed lines are for v\lv2 = 5, dotted lines are for v\/v2 = 0.2. In each panel, the upper family of curves are for 
1^20 = 10^*, middle for U2o = 10^^, and the fastest decaying curves are for i/20 = 10~^. The time units are such that |f2p|^^ = 1 at /? = 2. 
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Fig. 4. — Precession rate at i? = 2 for the models with S = 3/2. Lines correspond to the same models as for Fig. 3. The long dashed line 
shows the predicted Lense-Thirring precession rate at this radius. Note that the time axis is reduced as compared to Fig. 3. 
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Fig. 5. — Decay of the warp calculated using the linearized twist equations (solid lines) compared to the numerical solution of the full 
equations (dashed lines). The upper panel is for <5 = 3/2, j^io = 1^20 = 10~^, the lower panel is for the same parameters except that 
1^10 = 1^20 = 10"^. 



12 



Lense-Thirring precession of accretion disks 



1 1 1 1 r 



0.1 E 



10- 



10- 



10 



-4 



10 



-5 — - 



10 



10 



-6 



-7 _ 




10 



R 



Fig. 6. — Comparison of calculations using a Newtonian point mass potential (dashed lines) with a Paczynski-Wiita potential (solid lines). 
The parameters are as for Fig. 1. There are differences in the shape of the warp at late times, but the damping rates are almost identical. 



